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Detailed comparisons of the predictions of the Relativistic Form Factors (RFFs) and Modified Form Factors 
(MFFs) and their advantages and shortcomings in calculating elastic scattering cross sections can be found 
in the literature. However, the issues related to their implementation in the Monte Carlo (MC) sampling for 
coherently scattered photons is still under discussion. Secondly, the linear interpolation technique (LIT) is a 
popular method to draw the integrated values of squared RFFs/MFFs (i.e. A(Z, v?)) over squared momen- 
tum transfer (v? = v\, , Vg 9 ). In the current study, the role/issues of RFFs/MFFs and LIT in the MC 

sampling for the coherent scattering were analyzed. The results showed that the relative probability density 
curves sampled on the basis of MFFs are unable to reveal any extra scientific information as both the RFFs 
and MFFs produced the same MC sampled curves. Furthermore, no relationship was established between 
the multiple small peaks and irregular step shapes (i.e. statistical noise) in the PDFs and either RFFs or 
MFFs. In fact, the noise in the PDFs appeared due to the use of LIT. The density of the noise depends 
upon the interval length between two consecutive points in the input data table of A(Z, v?) and has no sci- 
entific background. The probability density function curves became smoother as the interval lengths were 
decreased. In conclusion, these statistical noises can be efficiently removed by introducing more data points 
in the A(Z, vf) data tables. 

Keywords: photons' elastic (Rayleigh) scattering; Monte Carlo techniques; probability distribution 
Function (PDF); relativistic form factors; modified form factors 



INTRODUCTION 

The application of Rayleigh (elastic) scattering by bound 
atomic electrons as being one of the fundamental modes of 
photon interaction at low energy X-rays and soft y-rays, to 
the different fields makes it of enhanced interest [1]. 
Furthermore, the availability of synchrotron sources and 
emerging development in the X-ray and y-ray regime has 
also boosted the demand for its accurate theoretical predic- 
tions [1, 2]. In the Rayleigh scattering, the whole atom acts as 
a primary target for the incoming photons while the electrons 
of the atom must be considered bound [3-5]. The target atom 
that has a relatively heavy atomic mass recoils to conserve 
momentum before and after scattering. Due to the effectively 



heavy atomic mass, the photon changes its direction only 
with negligible energy transfer to the target atom and retains 
the same energy after scattering [5]. Along with the atomic 
Rayleigh scattering by its bound electrons, nuclear Thomson, 
Delbriick and nuclear resonance scattering also contribute to 
the elastic scattering of photons by atoms [1]. 

Primarily, knowledge about Rayleigh scattering of photons 
by bound atomic electrons can be obtained using low energy 
X-rays and soft y-rays, in which the contribution from bound 
electrons of a neutral atom remains dominant. Now, if the 
electrons are considered free (i.e. bonding is neglected) the 
photon scattering will follow the famous Thomson scattering 
amplitude [6]. The differential Thomson cross-section for un- 
polarized photons is given by: 
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Here, r 0 = 2.82 x 10~ 13 cm is the classical electron radius [6, 7]. 
The differential cross-section for the coherent scattering may 
be obtained by multiplying the Thomson amplitudes and the 
square of the scattering atomic form factor F(Z, E, 0). The dif- 
ferential cross-section of coherent scattering within the atomic 
form factor approximation is given by: 



do- coh (Z,E, 6) dcr th {6) 



dO 



F 2 (Z,E,v) 



(2) 



Here, Z is the atomic number of the scattering element, E is 
the energy of the incident photon and v = A. _1 sin 9/2 (i.e. X 
(A) is the incident photon wavelength), having units of inverse 
length. The F(Z, E, v) can be written in terms of Relativistic 
Form Factors (RFFs) or Modified Form Factors (MFFs), and 
Anomalous Scattering Factors (ASFs) corrections as [3, 8]: 

F 2 (Z,E,v) = {F 0 (Z,v)+/'(Z,E)) 2 +/" 2 (Z,E) (3) 

Here, Fo(Z, v) represents the RFF/MFF, which is the 
Fourier Transform of the atomic charge distribution that 
signifies the phase difference in scattering from the differ- 
ing positions of the charge. In general, it represents the 
probability that the Z electrons of an atom take up the 
recoil momentum without absorbing any energy. It is also 
accounts for the effect of interfering scattering amplitudes 
across the atoms [4], while/' and/" are the real and the im- 
aginary parts of the complex scattering amplitude of the 
ASFs, accounting for the photoelectric influence [6, 7, 9]. 
The angle-independent ASFs are assumed to be isotropic 
and in addition, they smoothly approach zero at 1.0 MeV 
and, therefore, its effect can be neglected [10, 11]. 
Neglecting /'(Z, E) and/"(Z, E) in Equation (3) is consid- 
ered a further acceptable approximation for the Monte 
Carlo (MC) sampling of angular distribution of coherently 
scattered photons. By using Equation (1), Equation (2) can 
be rewritten as: 

djc ° h ^ a ' 9) = | (1 + cos^sinO F 2 (Z, v) (4) 

The form factor, FF, (i.e. RFF/MMF) approximation is 
valid for photon energies much greater than the electron- 
binding energies. Considerable effort has been made to 
obtain the precise values of RFFs/MFFs. Extensive tabula- 
tions of F 0 (i.e. RFFs and MFFs) for all neutral atoms have 
been published in the last quarter of the 20th century [2, 5, 
6, 9, 12]. These are traditionally tabulated as a function of 
variable v (A -1 ). The FF approximation is better justified at 
the energies well above the K-shell photo-effect threshold 
and generally gives good predictions for the total cross- 
sections [2, 3]. The FF approximation also works well to 



predict small angle differential cross-sections for the photon 
energies well above the K-shell photo-effect threshold. On 
the other hand, the effects related to atomic structure (i.e. 
virtual excitation and ionization of atomic electrons, at com- 
parable energies to the atomic binding energies especially 
for the photon energies close to the edges) are neglected in 
the FF approximation. So it fails to predict any of these 
structures. Owing to the dominant contributions from virtual 
ionizations and excitations of electrons in the subshells near 
the edges, it performs much more poorly than for that of the 
other energies [2]. 

Franz introduced for the first time in 1936 an 
electron-binding correction to the atomic form factors and 
the resulting corrected form factors are commonly known 
as MFFs [1, 2, 6]. The MFFs have not been calculated dir- 
ectly from the total electron charge distribution like RFFs. 
The contribution from the electrons of each subshell of an 
atom have been calculated and summed. The better high- 
energy results for heavy elements have been produced from 
such calculations. Furthermore, by the introduction of 
electron-binding effects, the MFF approximation provides 
the correct relativistic high-energy limit, particularly in 
high-Z elements, for forward scattering. The approximation 
is also quite successful in reproducing Rayleigh amplitudes, 
even at relatively low energies, where the validity of the 
derivation is unclear [2, 3, 5, 6]. In general, the accuracy of 
the MFF increases with increasing photon energy and de- 
creasing binding energies of the electrons [5]. Like the 
RFF, the MFF approximation also breaks down badly for 
high energy photons and at large scattering angles. 
However, unlike RFF, the MFF approximation gives good 
predictions for all elements at photon energies well above 
the K-shell binding after getting scattered by zero angles 
[3]. A complete tabulation of the MFFs values based on 
self-consistent relativistic wave-functions, potentials (i.e. 
total-atom and K-shell MFF values are listed) and by con- 
sidering the spherical symmetric electronic charge distribu- 
tion of the complete subshells given by Schaupp D et ai, 
for all the elements in the Periodic Table [5, 12]. 

Using the complete table of RFFs/MFFs covering the 
total energy range is almost impossible, or at least cumber- 
some, especially if it is included in a large production code 
where the speed of execution is of the primary importance. 
From this point of view, it is valuable to have simpler and 
faster methods to evaluate these data tables. The most 
popular method in the literature so far is the linear interpol- 
ation technique (LIT) for drawing the integrated value of 
squared RFFs/MFFs over squared momentum transfer (i.e. 

A(Z, v?)forv 2 = v 2 , , v 2 9 ) using the available data 

tables in the code [13]. 

A lot of research and discussion can be found about the 
predictions and shortcomings of RFFs and MFFs [1, 2, 14], 
but here the impact of electron-binding energies on coher- 
ently scattered photons was studied using MC techniques 
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for the sampling of coherently scattered photons. The 
impact has been assessed through the application of MFFs 
as it has electron-binding corrections. Another reason to 
evaluate its impact is to produce the best fit of possible 
probability distribution functions (PDFs) to the real angle 
distribution function for the coherently scattered photons. 
As LIT is the most popular method of drawing the A(Z, v?) 
values for vf however, the number of data points and interval 
length between two consecutive data points of the A(Z,vf) 
can be an issue. Here, the factors related to the validity of 
LIT and the interval length between two consecutive data 
points of the A(Z, vf) have been reviewed. 



kdp/dv 2 = — \/(ka) 2 can be easily established. Now 
equation (8) can be rewritten as: 



/ 2\ / 2 "■O'cofAZi, a, p) 
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(Tcoh(Z, a) 
Equation (2) can be written as 

o- C oh(Z,a, p,)d^ = TTr^l + p 2 )^o(Z, v)d|x 
By combining equations (9) & (10) 



(9) 



(10) 



MATERIALS AND METHODS 
Monte Carlo simulation 

In the past decades, several MC codes based on the rejection 
technique have been developed. The codes are either designed 
for the coherent scattering or for the transport of photons, 
containing the coherent scattering as an interaction mode of 
the photons. MC techniques established previously are reused 
in this study. [4,7, 15]. 

According to these techniques, the probability of the 
photons scattered coherently into the polar angle interval 
d9 around 9 is given by: 



p(v 2 )dv 2 



p(6)d6 



da- coh (Z,E, 9)d6 

o~coh[E, 0) 
r 2 (1 + cos 2 6)F 2 (v, Z)2irsin9d0 
2 <r coh {E, 6) 



(5) 



Where, o coh is the total coherent cross section [7]. 

For MC purposes, let us introduce p. = cos 6 and 
a = E/niQC 2 , (i.e. rriQC 2 = electron rest mass energy). Hence, 
the PDF in terms of u, given in equation (5), can be written 
as: [4] 



p(p)dp 



da- coh (Z, a, p)dp 
(T C oh(Z, a) 



(6) 



By using the trigonometric relation (i.e. sin 2 (6/2) = 
1/2(1 — cosO)) and the energy & wavelength relation, the 
relation, v = sin (8/2) /A can be easily written in the follow- 
ing form: 



v 2 = (ka)-(l - p), 0 < v 2 < v 2 



(7) 



where k = m^c / h\[2 — 29.1445 cnT x and v 2 corre- 
sponds to the v max [13]. 

Equation (6) can be written as, 



p(v 2 )dv 2 =p m h(p) 
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dv 2 
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Now by introducing, A(Z,v 2 ) = J F 2 ) (Z 1 v)Z 2 dv 2 and 



A(Z, v 2 ) = JFq(Z, v)Z 2 dv 2 for arbitrary v (these values 
o 

can be easily obtained by using numerical integration)equa- 
tion (8) will become: 

, , 2Trr 2 0 z 2 A(Z, v 2 ) (\ + p 2 \ F 2 (Z,E,v)z-W 
p[y )dv — — — y 



If, 



(ka) cr coh (Z 7 a) \ 2 



2 ._ F 2 (Z,E,v) Z - 2 



A(Z,v 2 ) 



(12) 



e(v 2 )= ' 7;;' -X ,/(v 2 ) = (^f) and 

2^z 2 A(Z,v 2 ) A ( Z ' V ) V 2 / 

2 1 men equation (12) can be written as 

(ka) a coh (Z,a) 

p( V 2 )d V 2 = c Q f(v 2 )Q( V 2 )dv 2 (13) 

To assign v 2 with density Q(v 2 ), a random number % on (0 
to 1) can therefore be used as follows [4, 1\. 



€= 



Q(v 2 )dv 2 = A(Z,E,v 2 ) /A(Z,E 



The value of v being accepted with probability [4, 7]: 

, 9N (l + M 2 ) 
/(v 2 )=^^<l 



(14) 



(15) 



By using equation (7) the relations p = 1 — [v 2 / (ka) 2 ] 



The integrated values of squared MFFs (i.e. Amff(Z,v 2 )) 
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and squared RFFs (i.e. A RFF (Z,vf)) over v 2 for 
vf = vf, , V25/vg 3/ / 59 /v| 5 has been calculated for hydro- 
gen (H, Z=l) H-l, Germanium (Ge, Z = 32) Ge-32 and 
Lead (Pb, Z= 82) Pb-82 by numerical integration [4, 7]. The 
percentage differences in Amff(Z,v?) and Ar/t/?(Z,v?) for 
for H-l, Ge-32 and Pb-82 were calculated 



1; 



: V 



15 ' 



■) "59 



using the following relation: 



% difference = 



(A RFF (Z,vf) - A MFF (Z,vf)) 
{(A RFF (Z,vf) 
x 100 



A MFF {Z,vf))/2} 



(17) 



To find the impact of the electron-binding correction on the 
angular distribution of -scattered photons, first the tabulated 
values of Arfp\Z, vf) for vf =v\, , V59, were 



included in the simulation code of coherent scattering. In 
the second step, the tabulated values of A RFF (Z,vf) were 

replaced by A MFF (Z, vf) for vf = v\, , (i.e. 

those consider the electron-binding energies) [2] in the 
simulation code. In both cases the PDFs were constructed 
through the above-mentioned MC techniques. During 
the program execution, the required A(Z, vf) values and v 2 
values can be found out through LIT by using the available 
input data tables of the Am FF (Z , vf) / A RFF (Z , vf) for 

vf = Vj, , V2 5 /v53/ 59 /vg 5 inside the code [4, 7]. The 

validity/effects of LIT on the MC sampling for the coherent 
scattering and the interval length between two adjacent data 
points of Am FF (Z, vf) in its input data tables were also ana- 
lysed. In doing so, 25 (n = 25), 53 in = 53) and 85 (n = 85) 
data points of MFF for V from 0 to 16.0 A -1 were used to 

calculate A MFF (Z, vf) for vf = vf, , ^5/^3/^5. On 

the basis of these values, the PDFs have been constructed 
through our MC code. The study was performed for H-l, 
Ge-32 and Pb-82 to find out the relation of the 
electron-binding correction and the interval length effect to 
atomic number (i.e. Z). 



RESULTS AND DISCUSSION 

The calculated values of A MFF {Z, vf) and A RFF {Z, vf) (i.e. 
for vf = vf, , vf 9 ) for H-l, Ge-32 and Pb-82 are pre- 
sented in Fig. la-c. The figure shows that Am FF {Z, vf) 
have smaller values compared to A RFF (Z,vf) for Ge-32 
and Pb-82, while the reverse is true for H-l. However, in 
both cases, the difference appears to increase for a while 
and then become constant. The mean, standard deviation 
(SD), minimum and maximum were calculated in the per- 
centage differences of Aypp{Z, vf) and A RFF (Z,vf) for 
H-l, Ge-32 and Pb-82. The results have been summarized 
in Table 1. The mean values, -2.138 51 ± 0.924 88, 
1.050 46 ±0.515 04 and 2.005 74 ±0.193 71, have been 
reported for H-l, Ge-32 and Pb-82, respectively (see Table 1). 
Fig. Id shows the mean of the % differences between 



vf = 0, ,vf 9 and A RFF (Z, vf) for 

vf = 0, , vj 9 for H-l, Ge-32 and Pb-82. 

Figure 2a-f shows the sampled relative probability 
density based on AM FF {Z,vf) and A RFF (Z,vf) with 

vf = vf, , v 2 9 , for the photons, having 0.001 MeV 

and 0.005 MeV energies after getting scattered by H-l, 
Ge-32 and Pb-82. The inset of Fig. 2a-f shows the mag- 
nified view of the relative probability densities of the 
angular distribution of the scattered photons for a specific 
range of u. In the magnified views, the difference 
between the curves, constructed using RFFs and MFFs, 
can be observed for both energies. The difference is 
more visible for 0.001 MeV compared to 0.005 MeV and 
looks most prominent in case of H-l. Fig. 3a-f shows 
the relative probability density for the photons having 
0.01 MeV and 0.05 MeV energies after getting scattered 
by H-l, Ge-32 and Pb-82. Here, the difference between the 
scattered angle distribution using Amff\Z, vf) and 
A RFF (Z, vf) is only visible for H-l at 0.01 MeV. Similarly, 
Fig. 4a-f shows the relative probability density for the 
photons having 0.1 MeV and 0.5 MeV energies. There is 
no distinguishable difference between the scattered angular 
distribution using Am FF (Z, vf) versus A RFF [Z, vf) on the 
scale of this plot, as can be observed from the insets of this 
figure. In general, it has been found that the inclusion of 
electron-binding energies makes some effect on relative 
probability densities at low energy, and as the energy of the 
photon increases the difference decreases. At high energy 
(i.e. >1 MeV) they are indistinguishable. It has also been 
found that in case of light elements like H- 1 , the difference 
is more visible compared to high Z elements, i.e. Pb-82. 
However, the literature has considered the MFFs as an 
improved approximation of exact scattering amplitudes, as 
long as the MFFs and RFFs are of the same order of magni- 
tude and at large scattering angles [5]. Nevertheless, the 
results reveal that the difference between sampled angular 
distribution using AM FF \Z,vf) and A RFF (Z,vf) is very 
small in the energy range -0.001-0.5 MeV. No extra in- 
formation was obtained using MFFs instead of RFFs in 
the sampling of PDFs for coherent scattering in this 
energy range. 

In Figs 2-4, non-smooth curves with multiple small 
peaks and irregular/step shapes, more prominent for the 
low energies, have been produced. The relation of these 
multiple small peaks and irregular/step shapes to MFFs/ 
RFFs and Z were investigated from the results obtained. 
These shapes were present in the relative probability 
density curves constructed using both the AM FF {Z,vf) 
andA RFF (Z,vf) for all the elements under the study, and 
hence their interconnection has been ruled out. Next, its re- 
lation with the MC techniques used for the sampling of co- 
herently scattered photons was investigated. These may be 
appeared due to some shortcomings in the MC techniques 
and are considered as statistical noise in the PDFs. It is a 
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Fig. 1. Amff[Z, V?) and A RFF (Z, v?)for v\ =0, , v^ 59 calculated on bases of MFF and RFF respectively for (a) hydrogen 

(H, Z=l), (b) Germanium (Ge, Z = 32), (c) Lead (Pb, Z = 82), (d) mean difference in % of A MFF (Z,vj) and A RFF (Z, v?) hydrogen 
(H, Z = 1), Germanium (Ge, Z = 32) and Lead (Pb, Z = 82). 



well-known fact that a random number generator (RNG) is 
at the heart of MC techniques. However, no relationship 
was found between the noise and either the RNG or the 
number of the simulated photons. 

The next step was to find out the relationship of these 
multiple small peaks and irregular/step shapes with the 
LIT used inside the MC code. For this purpose, the rela- 
tive probability density was sampled for 25 (n = 25), 53 
(n = 53) and 85 (n = 85) data points of A j hff(Z, vj) for 
the V in the range -0-16.0 A -1 . The MC samplings were 
performed for photons of energy in the range -0.0025- 
0.01 MeV after getting scattered by H-l, Ge-32 and 
Pb-82. Figure 5a-d shows the relative probability density 
for the photons, having 0.0025 MeV, 0.005 MeV, 0.0075 
MeV and 0.01 MeV energies by using n = 25, n = 53 and 
77 = 85 for H-l. The difference can be seen between the 
curves, constructed by using A MFF (Z,vf^ based on n = 25, 
n = 53 and n = 85 for all the energies in the study. The dif- 
ference is most prominent for n = 25 having some irregular 



Table 1. Mean, Standard Deviation (SD), Minimum and 
Maximum of percentage difference of A MFF (Z,vf) and A RFF 

(Z,vf) for vf = 0, , i'59 calculated on bases MFF and 

RFF respectively. 



Element 


Mean ± SD 


Minimum 


Maximum 


H-l 


-2.138 51 ±0.924 88 


-3.533 23 


0.004 22 


Ge-32 


1.050 46±0.515 04 


0.530 29 


2.506 37 


Pb-82 


2.005 74 ±0.193 71 


1.583 21 


2.432 69 



shapes. These irregular shapes at « = 25 are converted to 
the step shapes for n = 53 while the density of these 
shapes was further reduced for n = 85 at 0.0025 MeV. 
Figure 5a-d shows that as the photon energy was 
increased from 0.0025 MeV to 0.01 MeV, the density of 
the step shapes apparently decreased, but in fact the de- 
crease was not related to the energy but to the difference 
in the plot scale for different energies. Meaningful results 
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Fig. 2. Angle distribution of coherently scattered photon for source of 1 million photons. Inset shows the magnified view of relative 
probability density in the specific u (i.e. cosG) range, (a) Scatter element = H, energy = 0.001 MeV, (b) Scatter element = H, 
energy = 0.005 MeV, (c) Scatter element = Ge, energy = 0.001 MeV, (d) Scatter element = Ge, energy = 0.005 MeV, (e) Scatter 
element = Pb, energy = 0.001 MeV, (f ) Scatter element = Pb, energy = 0.005 MeV. 



196 



W. Muhammad and S.H. Lee 



(a) 



020 



10 



0.16 



3 0.12. 



O 



• 0.08 



O 

a. 



0.04 




(c) 



0 054 



0.045 



,0.036- 



I 



,0 027 



^0 013 



0 009 




# «FF' 

* UFF 



Ge-32 



Photon Energy: 0 01 MeV 



0.000 




I ■ I ■ I ■ I ■ I ' I ■ I ■ I ■ I 

-1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 
Mi M.I] 



(o) 



0 040 



£•0.012 

"35 
c 

& 



0 024 



O 0 016 

K 



0 008 



0 000 




Pb-82 



• M t» OK c»5 



Photon Energy: 0.01 MeV 



4» 

* utr 



9 

* 
* 

* 




I • I • I • I ' I • I » I • 1 • I 1 

•1.0 -0 8 -0.6 -0 4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 

M. M.1] 



(b) 



0 8 



in 

4 
o 

rx 0.4- 



o 
or 



0 2 



0.0. 



-1.0 



(d) 



0.30 



0.27- 
0.24 
0.21 
| 0.18 
S 0.16 
O 0 12 -J 



I 0.09 

j 

: 0.06 
0 03 -| 
0.00 
-0.03 



-1.0 



(f) 



0 20 ■ 



>> 0.15 

'35 



Q 

n 

o 



0.10 



Q) 0.05 

K 



000 



• 1 0 



ftM IN 



» un 



H-l 



Photon Energy: 0.05 MeV 



■0 5 



0.0 

Mi M.1] 



— i — 

0 5 



1 0 




4> KFF 
* Uff 



Ge-32 



Photon Energy: 0.05 MeV 




— i — 

•0.5 



1 

0.0 

v. M.1] 



i 

0.5 



1.0 



9 to on <M 



Pb-82 



Photon Energy: 0.05 MeV • 




•OS 



0 0 

M. M.1] 



0 5 



1 0 



Fig. 3. Angle distribution of coherently scattered photon for source of 1 million photons. Inset shows the magnified view of relative 
probability density in the specific u (i.e. cos0) range, (a) Scatter element = H, energy = 0.01 MeV, (b) Scatter element = H, energy = 0.05 
MeV, (c) Scatter element = Ge, energy = 0.01 MeV, (d) Scatter element = Ge, energy = 0.05 MeV, (e) Scatter element = Pb, energy = 0.01 
MeV, (f ) Scatter element = Pb, energy = 0.05 MeV. 
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Fig. 4. Angle distribution of coherently scattered photon for source of 1 million photons. Inset shows the magnified view of relative 
probability density in the specific u (i.e. cos0) range, (a) Scatter element = H, energy = 0.1 MeV, (b) Scatter element = H, energy = 0.5 
MeV, (c) Scatter element = Ge, energy = 0. 1 MeV, (d) Scatter element = Ge, energy = 0.5 MeV, (e) Scatter element = Pb, energy = 0. 1 
MeV, (f ) Scatter element = Pb, energy = 0.5 MeV. 
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Fig. 5. Angle distribution of coherently scattered photon for source of 1 million photons by using 25 (n = 25), 53 (n = 53) and 85 
(n = 85) data points of AFF for V from 0 to 16.0 A" 1 for H as a scatter element, (a) energy = 0.0025 MeV, (b) energy = 0.005 MeV, 
(c) energy = 0.0075 MeV, (d) energy = 0.01 MeV. 



can be obtained by inter-compression of the curves 
sampled on the basis of n = 25, n = 53 and n = 85 at the 
same energy. As a result, a clear difference can be seen 
between the curves sampled on the basis of n = 25, n = 53 
and « = 85. Similarly, Fig. 6a-d shows the probability 
densities for scattered photons, having 0.0025 MeV, 
0.005 MeV, 0.0075 MeV and 0.01 MeV energies for 
Ge-32 as a scattering element, while Fig. 7a-d shows the 
probability densities for the scattered photons, having 
0.005 MeV, 0.0075 MeV, 0.01 MeV and 0.05 MeV ener- 
gies for Pb-82 as a scattering element by using 
Amff(Z, v?) based on n = 25, n = 53 and « = 85 respect- 
ively. Like Fig. 5a-d, the same pattern has been observed 
but the difference between n = 25, n = 53 and n = 85 was 
decreased as Z for the scattering element was increased 



(i.e. Ge-32 and Pb-82) in the same range of energy. Now, 
it is clear from these figures that these multiple small 
peaks and irregular/step shapes are related to LIT. As the 
interval length between the two consecutive data points of 
Amff(Z, v?) was reduced, the density of multiple small 
peaks and irregular/step shapes in the probability density 
function was also reduced. Finally, it can be concluded 
that the multiple small peaks and irregular/step shapes 
depend upon the interval length between the two consecu- 
tive data points of Amff(Z, vf) in the input data tables. 
As the interval length increases, the shapes density in 
the relative probability density curves for the scattered 
photons increases, and the curve becomes smoother as the 
interval lengths are decreased. A smooth curve is 
expected if a large number of data points with the 
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Fig. 6. Angle distribution of coherently scattered photon for source of 1 million photons by using 25 (n = 25), 53 (n = 53) and 85 
(n = 85) data points of AFF for V from 0 to 16.0 A -1 for Ge as a scatter element, (a) energy = 0.0025 MeV, (b) energy = 0.005 MeV, (c) 
energy = 0.0075 MeV, (d) energy = 0.01 MeV. 



smallest possible interval lengths are included in the 
code. Furthermore, these can be considered as a statistical 
noise in the PDFs because of their relationship to sam- 
pling techniques rather than scientific background. 

CONCLUSION 

In conclusion, the MFFs had some effect, more visible for 
low Z than high Z elements, on relative probability dens- 
ities at low energy. However, the effect is very small and 
both RFFs and MFFs produced the same MC sampled 
angular distribution results for coherently scattered photons. 
No extra information was revealed by using MFFs instead 
of RFFs. Both can be used for the general purpose MC 
simulation of coherently scattered photons. Furthermore, no 
relationship was established between the multiple small 
peaks and irregular/step shapes in the PDFs either with 
RFFs or MFFs, as a result their scientific background has 



been ruled out. In fact, these appeared to be due to the use 
of LIT inside the MC sampling and can be treated as statis- 
tical noise in the PDFs. The interval length between two 
consecutive points in the input data tables of A(Z,vf) 
increases as the density of the noise increases and vice- 
versa. Furthermore, it can be concluded that a smooth 
curve can be produced if a large number of data points 
with the smallest possible interval length are used in the 
code. Therefore, to remove the statistical noise, the intro- 
duction of more data points with smaller interval lengths 
between two consecutive points in the input data tables of 
A(Z, v?) can be an efficient tool. 
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